Fermionic superfluidity: From high Tc superconductors to ultracold Fermi gases 
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We present a pairing fluctuation theory which self-consistently incorporates finite momentum pair excitations 
in the context of BCS-Bose-Einstein condensation (BEC) crossover, and we apply this theory to high T c 
superconductors and ultracold Fermi gases. There are strong similarities between Fermi gases in the unitary 
regime and high Tc superconductors. Here we address key issues of common interest, especially the pseudogap. 
In the Fermi gases we summarize recent experiments including various phase diagrams (with and without 
population imbalance), as well as evidence for a pseudogap in thermodynamic and other experiments. 
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Superfluidity in ultracold trapped fermionic gases has been 
one of the most exciting, rapidly evolving research subjects in 
the past few years. UTS H B HI 01. It has captured the at- 
tention of both condensed matter and atomic physics commu- 
nity, as well as other fields such as nuclear and quark matter 
|7|]. Atomic Fermi gases have strong similarities to supercon- 
ductors. Breakthroughs in cooling techniques in recent years 
have made it possible to cool an atomic Fermi gas to quantum 
degeneracy temperature and thus form a superfluid. What is 
remarkable about the atomic Fermi gases is their tunability 
and controllability. Using a Feshbach resonance, one can tune 
the attractive two-body interaction from weak to strong, and 
thereby make a smooth crossover from a BCS type of super- 
fluid to a Bose-Einstein condensation (BEC) BCS-BEC 
crossover has been of interest to the condensed matter com- 
munity for a long time, and has been argued [9] to be relevant 
to high T c superconductivity. 

Most of the theories that are proposed to explain high T c 
superconductivity are based on the pairing of electrons, as in 
BCS theory. While the pairing mechanism in the cuprate su- 
perconductors is far from being clear, there is no ambiguity 
about the pairing interaction in atomic Fermi gases. The at- 
tractive pairing potential between the fermionic atoms of dif- 
ferent spins has been well characterized and can be computed 
from first principle quantum mechanics calculations. Here we 
presume that cuprate superconductivity originates from elec- 
tron pairing, so that studying superfluidity in atomic Fermi 
gases may provide insights into the high T c superconductivity 
or vice versa fToL fTlfl . 

Importantly, preformed pairs are an inevitable consequence 
as the pairing interaction increases. At the unitary scatter- 
ing limit, the superfluid transition temperature T c is very high 
in units of the noninteracting Fermi energy Ep, The atomic 
Fermi gases often have a relatively higher transition tempera- 
ture than the cuprates. For the latter T c < Q.\Ep is somewhat 
lower, owing in part to their low dimensionality. Except in the 
extreme BCS limit, preformed pairs contribute a fermionic ex- 
citation gap, i.e., a pseudogap, even in the normal state. Such 
a pseudogap is common to both high T c superconductors and 
atomic Fermi gases. 

In the absence of a population imbalance, superfluidity is 
associated with a gas of fermions with dispersion £k = (v — 
fi, subjected to an attractive interaction between two different 



spin states, C/k.k' = Uip^ip^i, where U < is the coupling 
strength. For the cuprates, = 2fn(2 — cosk x — cosk y ), 
and for atomic Fermi gases ek = fc 2 /2m; \l is the fermionic 
chemical potential, t\\ is the in-plane hopping integral, and we 
have set the lattice constant to 1. Here and throughout, we 
assume natural units and set H = fcs = c = 1. The function 
</?k reflects the pairing symmetry. For the cuprates, ip-^ = 
cos k x — cos k y ; for the Fermi gases, the interaction is short- 
ranged, and taken to be a contact potential so that ipk = 1. 
The Hamiltonian is given by 

H - fiN = ^£k a L a k<r 

k<r 

+ ^2 C/k . k ' a k+q/2T a -k+q/2i a -k'+q/21 a k'+q/2T' 
kk'q 

(1) 

where a^(a) is the creation (annihilation) operator of the 
fermions. This picture is often referred to as a one-channel 
model. The physics of a Feshbach resonance can be described 
by a two-channel model, which includes an open-channel, as 
in the one-channel model, and a closed-channel, which rep- 
resents the effects of di-atomic molecules of a different total 
spin. The two-channel Hamiltonian is given by 

k,(T q 

+ X! f/ ( k ' k ') a q/2+k,T a q/2-k,i a q/2-kM a q/2+k',T 
q,k,k' 

+ ]T (.9(k)6 q a q/2 _ k ,i a q/2+k,T + hx ') ' < 2 > 

q,k 

where (b) is the creation (annihilation) operator for the 
closed-channel molecules, whose dispersion is given by e q = 
q 2 /2M with M = 2m. Here v represents the magnetic de- 
tuning, which is used to tune the relative energy level splitting 
between the open- and closed-channels and thus change the 
overall effective interaction; the interaction increases as the 
field decreases. The last line in Eq. (0 represents the cou- 
pling between the two channels, with g(k) = gv?k- 

Current studies of atomic Fermi gases have concentrated 
on 6 Li and 40 K. For both gases, the Feshbach resonances 
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which are used to enhance the pairing interaction have a very 
big width. As a consequence, the closed-channel fraction 
becomes negligible 0, Q even though the closed channel 
makes a big contribution to the effective pairing interaction. 
For example, for 6 Li, both experiments J6[] and our theoretical 
calculations lfl2ll have demonstrated that the closed-channel 
population is extremely small in the unitary regime (of the or- 
der of 10~ 5 — 10~ 4 ) and remains small in the entire range of 
magnetic field (i.e., pairing strength) which is accessed exper- 
imentally. 

The net effect of this tiny closed-channel fraction is that, for 
many purposes, one can use the one-channel model, Eq. (Q]), 
provided one replaces U with the overall effective interaction, 
U e ff(Q) = U + g 2 D Q (Q),whemD {Q) = l/[£fi„ - e™ 6 - 
v + 2/i] is the bare propagator of the closed-channel bosons. 
Experimentally, the pairing interaction for atomic Fermi gases 
is usually parametrized by the dimensionless product, l/kpa, 
where a is the s-wave inter-fermion scattering length, and kp 
is the Fermi wavevector defined in the noninteracting limit at 
T = 0, satisfying E F = k%/2m. 

In this way both high T c superconductors and atomic Fermi 
gases can be described by the same Hamiltonian but with dif- 
ferent pairing symmetries and free-particle dispersions. The 
existence of a big pseudogap and the short coherence length 
in the cuprates suggest that they are in the intermediate regime 
between BCS and BEC [9]. This is often referred to as the 
"strongly interacting regime" in atomic Fermi gases, and is 
the most complex and interesting regime as well. 

Our pairing fluctuation theory uses a T-matrix formal- 
ism. Instead of provi ding a detailed derivation, which can 
be found in Refs. ITol ll ill , here we recapitulate the key in- 
gredients and observations of this formalism. Pairs consist 
of fermions which glue together via the attractive interaction, 
and the pair propagator represents a summation of multiple 
particle-particle scattering processes. Fermions acquire self- 
energy via creation and destruction of a pair and a hole. In 
BCS theory, such pairs would be just the Cooper pairs in 
the condensate, and they contribute a self-energy H SC (K) = 
— A 2 c Gq(— K)(Pk' wnere A sc is the order parameter and 
Gq(K) is the bare fermion Green's function. [Throughout, 
we use a four-vector notation as in Ref. [10], K = (k, iu> n ), 
J2k = Sk' etc -' anc * w n is the Matsubara frequency.] 

In general, finite momentum pairs contribute to the fermion 
self-energy via E pg (K) = J2q tpg(Q)Go(Q - K)ip^, where 
t pg (Q) = 1/[1 + U ef f(Q)x(Q)} is the T-matrix (pair propa- 
gator), with x(Q) = J2k G o{Q ~ K)G(K) being the pair 
susceptibility and G{K) the full fermion Green's function. 
Below T c , t pg (Q) diverges at Q = so that we can approxi- 
mate Go(Q — K) k, Go(—K) in S ps . As a result, we obtain 
X pg (K) « -A^Goi-K)^, with 

a^-Ew^ (3) 

Q 

Here the key observation is that, under a reasonable approxi- 
mation, the self-energy of the finite momentum pairs has the 
same form as 'E SC (K). Since these finite momentum pairs do 
not have phase coherence, their contributions lead to a pseu- 
dogap in the fermion dispersion. The total excitation gap A is 



given by 

A 2 = A 2 C + A 2 ff , (4) 

and the fermionic quasiparticle dispersion becomes = 
\/£k + A 2 ^ 2 .. The fact that the gaps add in quadrature can 
be understood from the observation that A 2 C and A 2 g are pro- 
portional to the density of the condensed and noncondensed 
pairs, respectively. 

The effective dispersion 51 q and chemical potential n pa i r 
can be determined by a Taylor expansion of t~g, which, after 
analytical continuation (if2 n — > + i0 + ), can be written as 

tpg 1 = Z(tt - £l q + fip a i r ), (5) 

where Z is the inverse residue [11], and we have neglected the 
small imaginary part [11]. The pairing instability condition, or 
equivalently the BEC condition for pairs, requires 

*m(0) = 0, (6) 

so that we have fi pa i r — below T c . One can readily see 
that A pg — > as T — > 0. Therefore, this finite T theory is 
consistent with the BCS-Leggett ground state [8], which has 
been widely assumed in the literature. 

The presence of the pseudogap necessarily suppresses the 
superfluid transition T c . This can be seen from the relation- 
ship between A and T, as determined by a mean-field BCS 
gap equation: T decreases as A increases, and T c simply cor- 
responds to the lowest T at which A = A pg . 

Our theory can be summarized in three self-consistent 
equations: the gap equation (O, the pseudogap equation (0), 
and the number equation 

n = 2^G{K). (7) 

K 

They can be used to solve for A, A p9 , A sc , // and T c as a 
function of U (or U e t /) or 1/ kpa and T. 

The present theory naturally leads to the following impor- 
tant observations, (i) The fundamental statistical entities in 
these superfluids are fermions. We measure the "bosonic" (or 
pair) degrees of freedom indirectly via the fermionic gap pa- 
rameter A(T). (ii) As we go from BCS to BEC, pairs will 
begin to form at a temperature T* above T c . This pair for- 
mation is associated with a normal state pseudogap. (iii) In 
general there will be two types of excitations in these BCS- 
BEC crossover systems. Importantly in the intermediate case 
the excitations consist of a mix of both fermions and bosons. 

It should be pointed out that effects of pairing fluctua- 
tions were first considered by Nozieres and Schmitt-Rink II 1 311 . 
However, they only included them in the particle number 
equation; the absence of these effects in the gap equations 
necessarily leads to internal inconsistencies, say, in computing 
the superfluid density. In contrast to the present approach, the 
NSR scheme also predicts that T c approaches its BEC asymp- 
tote from above for an s-wave superfluid in three-dimensional 
(3D) continuum space. The NSR treatment, due to its rela- 
tive simplicity, has been used by others in the recent literature 

dm. 
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Figure 1: Cuprate phase diagram showing A(0), T c , A pg (T c ), and 
T*, calculated for -U/4t = 0.047, and tx/t\\ = 0.003. Shown 
as symbols are experimental data. The normal, pseudogap, and su- 
perconducting phases are labeled with "Normal", "P.G.", and "S.C.", 
respectively. FromRef. 11611 . 



A self-consistent treatment of finite-momentum pair excita- 
tions is crucial. This is more so in atomic Fermi gases, where 
one cannot directly measure the temperature except in the ex- 
treme BCS and BEC limits. Quantitatively reasonable knowl- 
edge of temperature effects is necessary in order to determine 
temperature and related physical properties in the BCS-BEC 
crossover regime. 

We first apply our pairing fluctuation theory to high T c su- 
perconductors. In Fig. Q] we plot the calculated phase dia- 
gram for E^Si^CaC^Os+a (Bi2212) and YB^CuaOy-a 
(YBCO) single crystals, as well as present a comparison with 
experimental data. The horizontal axis is the hole doping con- 
centration per unit cell, x = 1 — n. Here we show theoretical 
calculations and experimental data for zero T excitation gap 
A(0), superconducting transition temperature T c , and pairing 
onset temperature T*. We also show the theoretical result for 
A(T C ), the pseudogap at T c . 

To obtain this phase diagram, we fit with one parameter 
—U/Ato such that it gives the right value of T c at optimal 
doping around x = 0.15. Here to is the bare in-plane elec- 
tron hopping matrix element, in the absence of the Coulomb 
repulsion. It is related to tn via tu rj tox. One would ob- 
tain essentially the same phase diagram, if one parametrized 
the variable attractive interaction by simply fitting T* (x) to 
experiment. Thus, the underdoped regime sees effectively a 
stronger pairing interaction, — U/4tn, It is evident that our 
simple model gives a (semi-)quantitatively satisfying phase 
diagram, in agreement with experiment. 

It should be noted that, in contrast to the s-wave, 3D Fermi 
gas, the tf-wave symmetry and the lattice periodicity in the 
cuprates makes T c vanish before the BEC regime can be 
reached as the couping strength increases 11711 . This intro- 
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Figure 2: Comparison of superfluid density between theory and ex- 
periment. Top panel: Theoretical calculation for underdoped cuprate 
superconductors; Bottom panel: Experimental data from Ref. 11811 in 
the extreme underdoped regime of YBCO. 



duces a lower critical doping x c in the phase diagram shown 
in Fig.Q] 

Other applications of BCS-BEC crossover theory to the 
cuprate superconductors include studies ifioll of the quasi- 
universal behavior of the normalized superfluid density 
n s (T) /n s (0) as a function of T/T c for various doping con- 
centrations. Experimental data from Ref. Ill 811 are plotted in 
the lower panel of Fig. [2] for a series of very underdoped 
cuprates. The universality apparent in these data presents 
a puzzle for conventional theories of the penetration depth 
which invoke only fermionic excitations of the condensate. 
This is because the fermionic contributions are expected to 
reflect the strong variations in the excitation gap A(x) asso- 
ciated with the different hole concentrations x, thereby lead- 
ing to highly non-universal behavior. In the theory plot of 
the upper panel (Figure [2]) bosonic or pair excitations con- 
tribute an extra mechanism for destroying superfluidity. They 
compensate for the decrease of the quasiparticle contribution 
(oc T c (x)/A(x)) as the system evolves from overdoping to 
underdoping. As a result, the slope of the superfluid density 
remains nearly doping independent, with only a small system- 
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Figure 3: (color) Phase diagram of 40 K. A contour plot of the mea- 
sured condensate fraction No/N as a function of l/k%a and effec- 
tive temperature (T/Tp)° is compared with theoretically calculated 
contour lines at N s /N = (black curve) and 0.01 (red curve). The 
overall trend of the experimental contour of Nq/N = 0.01 and the 
theoretical line for N s /N = 0.01 are in good agreement. The dashed 
line represents the naive BCS result, T c jT% 



Ref. 



repre: 



0.615e 7r/2fc K From 



atic variation with x. The comparison between theory and re- 
cent experiments from the UBC group 111 811 shows quite good 
agreement. Recently, the bosonic nature of the superfluid den- 
sity in underdoped cuprates has also been recognized by other 
groups dEi. 

We turn now to experiments in cold atom systems. Here 
we ifTTl l22tl need to include the trap potential V(r). To this 
end, we use a local density approximation (LDA) by replac- 
ing /i — > fi(r) = fi — V(r). Here /i = /i(0) is the global 
chemical potential and V(r) = muj 2 r 2 /2 for a harmonic trap 
with mean angular frequency uj. The number constraint enters 
for the entire trap, N = J d 3 rn(r). Below T c , there exists a 
superfluid core within r < R sc - For r > R sc , ^L pa ir < 0, and 
the gap equation has to be modified to reflect this non-zero 

Impair • tpg (0) — Z[lp a i r . 

Before we compare with experiment, we have to know the 
temperature of the system. Experimentally, temperature is 
measured at the BCS [ 1] or BEC end ll23ll . before or after the 
experiment is performed in the strongly interacting regime. 
This BCS or BEC state is adiabatically connected to the state 
where experiments are done via a slow magnetic field sweep. 
By including self-consistently finite momentum pair excita- 
tions, our pairing fluctuation theory enables us to calculate the 
entropy S(T) at an arbitrary magnetic field. Then we can map 
the physical temperature T onto that measured in the BCS or 
BEC limit, T°, or vice versa Q. 

In Fig. [3] we show the phase diagram for 40 K and compare 
our calculations with experimental data [2 if]. Here all temper- 
atures are measured in (or converted into) the free Fermi gas 
limit. In addition, Tp = k F /2m is the Fermi temperature at 
the trap center in the noninteracting limit. The experimental 
superfluid phase boundary is given by the 1 % condensate frac- 
tion contour line; it is hard to locate experimentally where the 
condensate vanishes precisely. This should be compared with 
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Figure 4: (color) Energy per particle as a function of T for 6 Li at 
unitarity. The fact that the experimental data (symbols) (and the two 
theoretical curves) for noninteracting and unitary Fermi gases do not 
merge until higher T* > T c is consistent with the presence of a 
normal state pseudogap. From Ref. 12511 . 



the (red) theory curve. Here we 12111 identify the condensate 
fraction with the superfluid density. The agreement between 
experiment and theory is reasonably good. It should be noted 
that both experiment and theory show a maximum around uni- 
tarity, l/k%a = 0, which has also been observed in 6 Li phase 
diagrams jj]. 

We next turn to a detailed comparison of theory and exper- 
iment for thermodynamics. Figure [4] presents a plot of en- 
ergy E as a function of T comparing the unitary and non- 
interacting regimes for 6 Li. The solid curves are theoretical 
while the symbols are experimental data lr25tl . A realistic 
Gaussian trap with the experimentally given trap depth was 
used. There has been a recalibration of the experimental tem- 
perature scale i25ll in order to plot theory and experiment in 
the same figure. The latter was determined via Thomas-Fermi 
fits to the density profiles [22]. To arrive at the calibration, 
we applied the same fits to the theoretically produced density 
profiles. 

Good agreement between theory and experiment is appar- 
ent in Fig. |4] It should be emphasized that there is no fitting 
parameter in our calculations. The temperature dependence 
of E reflects primarily fermionic excitations at the edge of 
the trap 12411 . although there is a small bosonic contribution as 
well. Importantly one can see the effect of a pseudogap in the 
unitary case. That is, the temperature T* is visible from the 
plots as that at which the non-interacting and unitary curves 
merge. This corresponds roughly to T* « 2T C . Evidence of 
a pseudogap at unitarity has also been observed more directly 
in a radio-frequency excitation gap experiment ll23ll . 

Finally, we turn to the effects of population imbalance, 
which have recently been one of the hottest subjects in atomic 
Fermi gases. We consider both the homogeneous system 1 26], 
as well as trapped case, which can be found in Refs. B271 12811 . 
We define n = n-t + n\, Sn = — n^, and the polarization 
p = 8n/n. There will be two number equations associated 
with the two spin species. The polarized superfluid phase is 
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Figure 5: (color) Phase diagram in the p — 1/fejra plane, show- 
ing where intermediate temperature superfluidity associated with the 
Sarma phase (shaded region) exists. The (red) curve labeled by 
d 2 £l(T c )/dA 2 — is determined at corresponding T c . The line 
defined by da/dT c — is given by the turning points (p, l/fc^o) 
where 1/kpa reaches a local extremum as a function of T c . At T = 
0, the entire region between the T C MF = and the d 2 £l/dA 2 = 
lines is unstable against phase separation. However, a A stable po- 
larized, intermediate temperature superfluid phase exists for the (yel- 
low) shaded region. Homogeneous superfluids in the dotted region 
are unstable at any temperature. From Ref. 13111 . 



pears consistent with experiment. More details can be found 
in Ref. fig] . 

Finally, we show in Fig.|6]the phase diagram at unitarity for 
a population-imbalanced Fermi gas in a harmonic trap with 
angular frequency u>. Here 1/kpa = 0. Phase separation (la- 
beled PS) occupies the lower T portion of the phase diagram, 
where the gap A jumps abruptly to zero at some trap radius. 
At intermediate T, there is a (yellow-shaded) Sarma phase, 
where A vanishes continuously within the trap. It evolves 
into a (dotted region) pseudogap (PG) phase as the superfluid 
core vanishes at higher T. A normal (N) phase without pairing 
always exists at even higher T. This phase diagram appears 
to be consistent with current experiments. One can easily see 
that the physics associated with population imbalance is much 
richer than in the absence of imbalance. 

In conclusion, in this paper we have addressed common- 
alities particularly associated with the phase diagrams of ul- 
tracold trapped fermionic gases and high T c superconductors. 
These common features revolve around the scenario [Till 
that the cuprates are somewhere intermediate between BCS 
and BEC. Importantly, this scenario has been directly real- 
ized in trapped Fermi gases. Here one sees considerable evi- 
dence for pre-formed pairs and the related fermionic pseudo- 
gap which appear reminiscent of their cuprate counterparts. 

This work was supported by NSF PHY-0555325 and NSF- 
MRSEC Grant No. DMR-02 13745. 



conventionally called the "Sarma" phase [29], although here it 
is generalized away from the BCS limit originally considered 
by Sarma. Details are given in Ref. |26]. 

In Fig. |5J we present the zero T phase diagram for a popu- 
lation imbalanced homogeneous Fermi gas in the p — 1/fcpa 
plane. Here T C MF is the solution of T c in a strict mean-field 
treatment. Q, is the thermodynamic potential. A superfluid 
solution is stable only if d 2 fl/dA 2 > 0. From this figure, 
we see that at any p ^ 0, there is a minimum threshold for 
1 / k,Fa (determined by T^ IF = 0) in order to have a super- 
fluid solution. The point at which the solution becomes stable 
(determined by d 2 il/dA 2 = 0) occurs in BEC regime. The 
(blue) dashed line indicates where T c vanishes. The dotted re- 
gion (as well as the shaded region) in the figure indicates an 
unstable (Sarma) superfluid phase at T = 0. The stable phases 
in this region are presumed to be either phase separated (PS) 
or possibly a Fulde-Ferrell-Larkin-Ovchinnikov state li3Qfl . 

In addition to the zero T phases, we also show in Fig. [5] 
where the system is stable at T c . A stable finite tempera- 
ture polarized superfluid exists in the (yellow) shaded region. 
This is a very interesting prediction of our theory which ap- 
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Figure 6: (color) Phase diagram of a population-imbalanced Fermi 
gas in a harmonic trap at unitarity. The solid lines separate different 
phases. Above the (red) dashed line but within the PS phase, the su- 
perfluid core does not touch the domain wall. Here "PG" indicates 
the pseudogapped normal phase. The black dot labeled "TCP" indi- 
cates the tricritical point. 
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